DISCRETE FOURIER ANALYSIS ON A DODECAHEDRON AND 

A TETRAHEDRON 



HUIYUAN LI AND YUAN XU 



Abstract. A discrete Fourier analysis on the dodecahedron is studied, from 
which results on a tetrahedron is deduced by invariance. The results include 
Fourier analysis in trigonometric functions, interpolation and cubature formu- 
las on these domains. In particular, a trigonometric Lagrange interpolation 
on the tetrahedron is shown to satisfy an explicit compact formula and the 
Lebesgue constant of the interpolation is shown to be in the order of (logn) 3 . 



1. Introduction 

It is well known that Fourier analysis in several variables can be developed based 
on the periodicity defined by a lattice, which is a discrete subgroup defined by ^4Z d , 
where A is a nonsingular d x d matrix. A lattice L := AZ d is called a tiling lattice 
of M. d if there is bounded set fl that tiles M. d in the sense that f2 + L = R d . Let L 
be a tiling lattice and L 1 - := A~ tr Z d be its dual lattice; then a theorem of Fuglede 
[5] states that the family of exponentials \^ 2lxla x : a G L } forms an orthonormal 
basis for L 2 (fl). The Fourier expansion on fl is essentially the usual multivariate 
Fourier series under a change of variables x i— > A~ x x. 

One can also develop a discrete Fourier analysis associated with a lattice, starting 
with a discrete Fourier transform based on L 1 - , which has applications in areas such 
as signal processing and sampling theory (see, for example, [H HQ]). Recently 
in [5], we studied the discrete Fourier transform and used it to derive results on 
cubature and trigonometric interpolation on the domain O, both are important 
tools in numerical computation and approximation theory. The simplest domain for 
the tiling lattice is the regular hexagon, which has the invariance of the reflection 
group A2 • The fundamental domain of the hexagon under A2 is an equilateral 
triangle. A detailed study of the discrete Fourier analysis is carried out on the 
hexagon and on the triangle in [S]. The invariant and the anti-invariant projections 
of the basic exponential functions are analogues of cosine and sine functions on 
the triangle, which have been studied previously in [5j 111], Explicit and compact 
formulas are derived for several cubature formulas and interpolation functions in 
In particular, we found a compact formula for the Lagrange interpolation by 
trigonometric functions that interpolates at X n := {(^, ^) : < i < j < n} on the 
triangle T := {(x, y) : x, y > 0, x + y < 1} and proved that its Lebesgue constant 
is in the order of (logn) 2 . The result on interpolation is noteworthy since it is 
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in sharp contract to the algebraic polynomial interpolation on X n , which has an 
undesirable convergence behavior. 

The purpose of the present paper is to carry out a similar analysis on R 3 for a 
tetrahedron, also called a simplex in R 3 . For this we work with the face-centered 
cubic (fee) lattice, which has the symmetry of reflection group A3. The domain 
that tiles R 3 with the fee lattice is the rhombic dodecahedron (see Figure 3.2), 
whose fundamental domain under .A3 is a regular tetrahedron. We shall develop in 
detail a Fourier analysis on these two domains, study analogues of cosine and sine 
functions as in the case of hexagon, and establish compact formulas for discrete 
inner product, cubature formulas and Dirichlet kernels. 

Just as in the case of the regular hexagon [9l [11] , the analysis on the rhombic 
dodecahedron and the tetrahedron is carried out using homogeneous coordinates of 
R instead of in R 3 . This has the advantage that our formulas are more symmetric 
and the symmetry of the domain becomes more transparent. The Fourier transform 
on the dodecahedron as well as the generalized cosine and sine functions were 
studied earlier in [12] using a homogeneous coordinate system in R 6 . We choose 
our homogeneous coordinates in R 4 since A3 can be regarded as a permutation 
group on four elements. 

It should be pointed out that the development on a specific domain does not 
follow immediately from the general theory. To tile the space without overlap, the 
domain Q can only include part of its boundary. For the discrete Fourier analysis, 
this fact causes a loss of symmetry; for example, for the rhombic dodecahedron, 
the discrete Fourier transform is defined using only part of the boundary points. In 
order to obtain results that are symmetric, we have to modify definitions to include 
all boundary points, which can be delicate if the orthogonality is to be preserved. 
The difficulty lies in the congruent relations of the boundary. In order to understand 
the periodicity based on the rhombic dodecahedron, we need to understand the 
congruence of the boundary under translation by the lattice. Furthermore, in order 
to transform results from the rhombic dodecahedron to the tetrahedron, we need to 
understand the action of A3 on the boundary. The complication of the congruence 
of the boundary is also one of the main reasons why we restrict ourself to R 3 instead 
of dealing with lattices on R d that are invariant under Ad (see [2]) for all d > 3. 

One of our main results is a compact formula for the Lagrange interpolation 
based on the regular points on the tetrahedron, whose Lebesgue constant is shown 
to be in the order of (logn) 3 . Again, this is a result in sharp contrast to interpo- 
lation by algebraic polynomials on the same set of points. Interpolation by simple 
functions is an important tool in numerical analysis that has a variety of appli- 
cations. For interpolation on the point sets in several variables, little results are 
known if the point sets are not of tensor product type. Moreover, most studies 
consider mainly interpolation by algebraic polynomials, which face the problem of 
choosing interpolation points, as equally spaced points do not yield favorable re- 
sults [1] . Our study in [9] and in the present paper demonstrates that interpolation 
at equally spaced points on the triangle and on the tetrahedron can be solved with 
trigonometric functions: the interpolation can be carried out by compact formu- 
las that offer fast computation, and the convergence behavior is as good as can 
be expected since a Lebesgue constant of (log n) d for interpolation in R d is about 
optimal. In the present paper, we concentrate on theoretic framework that leads 
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to trigonometric interpolation on the tetrahedron, numerical study is left out for a 
future work. 

The paper is organized as follows. In Section 2 we sum up results from the 
general theory of discrete Fourier analysis associated with lattice. The analysis 
on the rhombic dodecahedron will be carried out in Section 3, including a detailed 
study on the congruence of the boundary. Results on the tetrahedron are developed 
in Section 4, including generalized sine and cosine functions. 

2. Discrete Fourier analysis with Lattice 

In this section we recall results on discrete Fourier analysis associated with lat- 
tice. Background and the content of the first subsection can be found in [21 [31 [71 [TU] . 
Results in the second subsection are developed in [5] . We shall be brief and refer 
the proof and discussions to the above mentioned references. 

2.1. Lattice and Fourier series. A lattice L of R d is a discrete subgroup that 
contains d linearly independent vectors, 

L := {fcidi + k 2 a 2 H + k d a d : k t G Z, i = 1, 2, • • • , d} , 

where eti, • • • , ad are linearly independent column vectors in M. d . Let A G M. dxd be 
the matrix whose columns are a\, ■ ■ ■ , ad- Then A is called a generator matrix of 
the lattice L. We can write L as La and a short notation for La is ^4Z d ; that is 

L A = AZ d = {Ak: fee Z d } . 

Throughout this paper, we shall treat a vector in the Euclidean space as a column 
vector whenever needed. As a result, x ■ y — x tr y, where x tr denotes the transpose 
of x. The dual lattice L of L is given by 

L 1 - := {x G M. d : x ■ y £ Z for all y E L) . 

where x ■ y denotes the usual Euclidean inner product of x and y. The generator 
matrix of L is A~ tr . 

A bounded set fi C K d is said to tile M. d with the lattice L if 

Xn(x + oi) = 1, for almost all x G M. d , 

where XQ denotes the characteristic function of f2, which we write as f2 + L = M. d . 
Tiling and Fourier analysis are closely related as demonstrated by the Fuglede 
theorem. Let f(x)dx denote the integration of the function / over f2. Let (•, -)q 
denote the inner product in L 2 (£l), 

(2.1) (f,9)a~^ff{x)g(x)dx, 

where |Q| denotes the measure of fl. The following fundamental result was proved 
by Fuglede in [5]. 

Theorem 2.1. Let 11 C M. d be a bounded domain and L be a lattice ofM. d . Then 
SI + L = M. d if and only if | e 27rlQ a; : a g L^ - } is an orthonormal basis with respect 
to the inner product (|2.1| . 



4 



HUIYUAN LI AND YUAN XU 



The orthonormal property is defined with respect to the normalized Lebesgue 
measure on fi. If L = La, then the measure of is equal to det(A tr A) . Fur- 
thermore, we can write a G L\ = A~ tr Z d as a = A~ tr k with k e Z d , so that 
a ■ x = k tr A~ l x. Hence the orthogonality in the theorem is 

(2.2) - 1 [ e 2mkt ' A ~ lx dx = S k0 , ke Z d . 

V ' ^/det(A^A) Jn 

The set f2 is called a spectral set (fundamental region) for the lattice L. If L = La 
we also write — £Ia- 

A function / G L 1 (Qa) can be expanded into a Fourier series 

fix) ~ Y c k e 2 ^ kt ' A "\ c k = - 1 / f(x)e-^ ik% ' A ~' x dx. 

The Fourier transform f of a function defined on L 1 (W i ) and its inversion are 
defined by 

/(£) = / f(x)e- 2 ™^dx, f{x) = [ /i(;«-' ; '*rf(. 

JR d JR d 

Our first result is the following sampling theorem (see, for example, [71 [TU]V 

Proposition 2.2. Let be the spectral set of the lattice AZ d . Assume that f is 
supported on fl and f G L 2 (Q). Then 

/(*)= E f(A- tr k)$ n {x-A- tr k) 



in L 2 (fl), where 



^/det(A^A) Jn 



This theorem is a consequence of the Poisson summation formula. We notice 
that 

$ Q (A- tr j) = 5 j, foralljeZ rf , 
by Theorem 12. l\ so that <£>n can be considered as a cardinal interpolation function. 

2.2. Discrete Fourier analysis and interpolation. A function / defined on R d 
is called periodic with respect to the lattice AZ d if 

f(x + Ak) = f(x) for all k G Z d . 

The spectral set f2 of the lattice AZ d is not unique. In order to carry out the 
discrete Fourier analysis with respect to the lattice, we shall fix f2 such that 
contains in its interior and we further require that f2 tiles M. d with La without 
overlapping and without gap. In other words, we require that 

(2.3) E Xn(x + Ak) = 1, for all .t 6 

feez d 

For example, we can take i7 = [— |, ^) d for the standard cubic lattice Z d . 

Definition 2.3. Let A and B be two nonsingular matrices in E dx<i ; VIa and VIb 
satisfy (|2.3|) . Assume all entries of N := B tr A are integers. Define 

A N := {k G Z d : B~ tr k G ft A } and := {k E Z d : A~ tr k G Q B } ■ 
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Two points x,y €M. d are said to be congruent with respect to the lattice AZ, d , 
if x — y G AZ d , and we write x = y (mod A). The following two theorems are the 
central results for the discrete Fourier transform. 



Theorem 2.4. Let A, B and N be as in Definition \2.3\ Then 

1 J2 -2-U«;v-S-_ f 1 ' */fc = (mod7V tr ), 



I det(JV) I *-f I 0, otherwise 



Theorem 2.5. Let A, B and N be as in Deftnition \2.3[ Define the discrete inner 
product 

/or /, g G £7(0,4 ), fAe space o/ continuous functions on Qa- Then 
(2-4) (/, .g}^ - (/, 5 ) 

for all f, g in the finite dimensional subspace 

H N := span {<fe : cj> k [x) = e 2 ™^ 1 *, k e A^} . 

Let |.E| denote the cardinality of the set E. Then the dimension of TCn is Ajy|. 

Let Xjv/ denote the Fourier expansion of / 6 (7(0^) in Hat with respect to the 
inner product (•, -}jv- Then, analogous to the sampling theorem in Proposition ^. 21 
Ijv/ satisfies the following formula 

where 



~ |det(AT)| ^ 
feeAt 



e 27rife"A 'i 



The following theorem shows that In} is an interpolation function. 

Theorem 2.6. Let A, B and N be as in Definition \2.3[ If A' N — Ajv», then I^f 
is the unique interpolation operator on N in Hn; that is 

l N f(B- tr j) = f{B-«j), Vj e A N . 

In particular, \An\ — |AL|. Furthermore, the fundamental interpolation function 
$q b satisfies 

k£Z d 

The above results have been used to develop a discrete Fourier analysis on a 
hexagon in [9]. In the following section, we apply it to the rhombic dodecahedron. 

3. Discrete Fourier analysis on the rhombic dodecahedron 

In this section we develop a discrete Fourier analysis on the rhombic dodecahe- 
dron. It contains five subsections. 
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Figure 3.1. Face-centered cubic (fee) lattice. 




3.1. Face-centered cubic lattice and Fourier analysis. We consider the face- 
centered cubic (fee) lattice given in Figure 3.1. 

Just like the hexagon lattice, the fee lattice offers the densest packing of R 3 with 
unit balls, which is the so-called Kepler's conjecture and proved recently in [6]. 
The generator matrix A of fee lattice is given by of 



.4 



The domain that tiles R with fee lattice is the rhombic dodecahedron (see Figure 
3.2). Thus, the spectral set of fee is 

fi = {x e R 3 : -1 < Xj ; ± x t < 1, 1 < i < j < 3}. 

The strict inequality in the definition of f2 reflects our requirement that the tiling 
of the spectral set has no overlapping. 

Motivated by the study of [HI El [H], we shall use homogeneous coordinates 
t = (tt,tz, t3, ti), where t\ + t% + t 3 + t± = 0, in R 4 for our analysis on the rhombic 
dodecahedron in R 3 . The advantage is that our formulas become more symmetric 
and the symmetry becomes more transparent under the homogeneous coordinates. 
Throughout the rest of this paper, we adopt the convention of using bold letters, 
such as t, to denote the points in the space 

R% := {t = (ti, ta, t 3 , U) E R 4 : h + t 2 + t 3 + t± = 0} . 

In other words, the bold letters such as t and k will always mean homogeneous 
coordinates. The transformation between x G R 3 and t £ is defined by 



(3.1) 




Xl = t 2 + t 3 
x 2 =h+ t 3 
x 3 =t 2 + t 1 
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Figure 3.2. Rhombic dodecahedron 



H 



and ti — — 1\ — t2 — t 3 . Let us denote by H and U the matrices 

(- 1 1 

1 -1 

1 1 

V-i -1 -1/ 

respectively. The columns of the matrix U are orthonormal and U tr U 
then have A = U tr H and the inverse transform of (13.111 is 








°\ 





1 











1 


V-i 


-1 


-v 



and U = — 
2 



1\ 
1 

-1 



(3.2) 



Ux 



I. We 



t\ = \{-Xi +X 2 + X 3 ) 
t-2 = 5(^1 - X 2 + X 3 ) 

h = \{x\ +x 2 ~ x 3 ) 
t± = \{~x l -x 2 - x 3 ). 

In the homogenous coordinates, the spectral set is fin := {t = Ux : x £ Vt] 
which, upon using (|3.2p . results to 

(3.3) VL H = {t e lg : -1 < U - tj < 1, 1 < i < j < 4} . 

Figure 3.3 shows again the rhombic dodecahedron with vertices labeled in the ho- 
mogeneous coordinates. 
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Figure 3.3. Rhombic dodecahedron labeled in homogeneous coordinates. 



Under the change of variables (|3.2[) , the integral over is given by 

(M ra/«. /(t) *-i5i/ /Wfc 

By Fuglede theorem, forms an orthonormal basis in L 2 (17). 

We would like to reformulate the exponential functions e 27 ™ fe A 1 so that they are 
indexed by homogeneous coordinates. For this purpose, we denote by 

Z^ := Z 4 n M.% = {k G Z 4 : fci + fc 2 + fc 3 + fc 4 = 0} 

the set of integers in homogeneous coordinates and introduce the notation 

(3.5) H := {k G Z% : k t = k 2 = k 3 = k 4 (mod 4)} . 

The definition of the matrix H shows that if k G Z 3 then Hk G 1i A H . For k G Z 3 , 
set j = 4H(A tr A)~ 1 k G Z^. A quick computation reveals that j G H. Moreover, 
given j G H, it follows from the fact A = U tr H that k = jA tr U tr j, which is easily 
seen to be in Z 3 . Furthermore, we have 

k^A^x = \fUAA~ x x = \fUx = ij • t. 

Consequently, we can index the exponentials by j G H and the exponent 2irik tr A~ 1 x 
becomes • t. Let us introduce the notation 

(3.6) 0j(t) :=e^\ j G EL 

Then, using (|3.2p and recalling (|2.2p . the Fuglede Theorem becomes the following: 
Proposition 3.1. For k, j G HI, 

(0k,</>j) = o / <Mt)<Aj(t)dt = 5 k ,j. 

2 JS2h 

Furthermore, {<p- } : j G H} is an orthonormal basis of L 2 (Q,h)- 

Given / defined on £1, the mapping (|3 . 2[) shows that /(x) = f(U ir t) = g(t) is the 
function in homogeneous coordinates. Since A = U tr H, a function / being periodic 
with respect to the lattice Alfi becomes, in homogeneous coordinates, the following 
definition: 
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Definition 3.2. A function f is H -periodic if it is periodic with respect to the fee 
lattice; that is, /(t) = f(t + Hk) for x G £Ih and k G Z 3 . 

Using the explicit form of the matrix H, it is easy to see that the following holds: 

Lemma 3.3. A function /(t) is H -periodic if and only if 

/(t) = /(s), t-seZ%, 

or equivalently 

/(t + ke itj ) = /(t), fceZ, l<Kj<4, 

where e,-j :— e; — ej and {ei, ei, e3, 64} is i/ie standard basis ofM 4 . 

Evidently, the functions </>j(t) in ()3.6|) are ii-periodic. Furthermore, Proposition 
13.11 shows that an i7-periodic function / can be expanded into a Fourier series 

(3.7) /~^/k<£k(t), where h ■= \ f /(t)^_ k (t)dt. 

3.2. Boundary of the rhombic dodecahedron. In order to carry out the dis- 
crete Fourier analysis on the rhombic dodecahedron, we need to have a detailed 
knowledge of the boundary of the polyhedral. 

We use the standard set theory notations dil, ft and fl to denote the boundary, 
the interior and the closure of f2, respectively. Clearly 57 = tt° U dil. A rhombic 
dodecahedron has 12 faces, 24 edges and 14 vertices. Since we will consider points 
on the boundary, we need to distinguish a face with its edges and without its edges, 
and an edge with its end points and without its end points. In the following, when 
we say a face or an edge, we mean the open set, that is, without its edges or end 
points, respectively. 

We shall work with homogeneous coordinates. To describe the boundary of Qh 
we set N 4 := {1,2,3,4}. For i, j G N 4 and i ^ j, define 

F itj = {t G U H : U - tj = 1}. 

There are a total 2(*) = 12 distinct F^j and it easy to see that each Fij stands for 
one face, with its edges, of the rhombic dodecahedron. 
For nonempty subsets /, J of N4, define 

Qij ■= P| F hj = {t G U H : tj =U-l, for all i £ I. j G J} . 

ieijeJ 

Lemma 3.4. Let I,J.Ii,Ji be nonempty subsets 0/N4. Then 

(i) fl^j = if and only if I fl J 7^ 0. 

(ii) fl IuJl n fl l2j2 = Q, ItJ ifh U I 2 = I and J 1 UJ 2 = J. 

Proof. It is obvious that Qij ^ if / n J = 0. On the other hand, if / n J ^ and 
i = j G / H J, then U — tj = 7^ 1, which shows that ili.j = 0. This proves (i). 

If either h n J\ / or I 2 n J 2 ^ 0, then fli uJl l~l fli 2 ,j 2 = fij.j = by (i). If 
I\ n Ji = I2 n J2 = and i v G I v , j v G J v for 1/ = 1, 2, then 

til — ij 2 = ^2 ~l~ 1 — 1 and — 1 < £j 2 = £j 2 — 1, 
which implies — ti 2 =0 so that tj 1 = tj 2 = t% x — 1 = U 2 — 1 and proves 

nj,,* nfti 2 ,j 2 = (i It j. " ' □ 



10 



HUIYUAN LI AND YUAN XU 



Edges are intersections of faces and vertices are intersections of edges. The 
Lemma [3.4l gives us information about the intersections. To make clear the structure 
of the boundary dfln, we introduce the notation 

K = {(I,J) :/, JcN 4 ; I n J = 0} , 

/C = {(/, J)GfC: i < j, for all (i, j) G (I, J)} . 

Definition 3.5. For (/, J) € /C, f/ie boundary element Bu of the dodecahedron, 

B hJ := {t G : t g nj,,^ /or a/Z J x ) € /C wiift |J| + |J| < |ii| + |Ji|}, 

is called a face if \I\ + \J\ = 2, an edge if \I\ + | J\ = 3 and a vertex if \I\ + | J| = 4. 

For the faces and the edges, the boundary elements represent the interiors. In 
fact, it is easy to see that = F°j and, for example, -B{i}.{j,fc} = {Fij ^Fi,k)° 

for distinct integers i, j, k £ N4. 

Furthermore, for < i, j < i + j < 4, we define 

AC ij 1 :={(/, J)e/C: |I|=i, I J|=j}, |J 

(3.8) ■ • , , 

/C J := {(J, J) G /Co : |/| = i, \J\ = j} , Bj' := (J 

Note that is the union of boundary points in those Bj j for which \I\ = i and 
\J\=J- 

Proposition 3.6. Let (I, J) G K and (h,J\) G /C. 

(i) Bjj n flj^jj =%,ifI^I 1 andJ^J l . 

(ii) n H \tt° H = U(/,j)eJc = Uo<i,i<i+j<4 fiI J • 

(iii) fa \ Q° H = U(/ lJ)e) Co = Uo ,, , , ; 



Proof. If Bi tJ nBi uJl ^ then |J| + |J| = |Ii| + |Ji|. Moreover, if t £ Bi tJ nBj uJl 
then t £ f^j n Oji.jj = Oju/i,juJi by Lemma T3. 41 which implies that |/| + | J\ > 
\I U 7i| + I J U Ji\. Thus we must have I = 7i = IU h and J = J x = J U Ji, which 
contradicts the assumption and proves (i). 

To prove (ii), we define for t £ <9£!#, I = {i G N4 : 3j £ N4 such that i; — tj = 1} 
and J = {j G N4 : 3i G N 4 such that U - tj = 1}. Clearly (I, J) G /C and t G P/,./, 
which proves the first equal sign of (ii). The second equal sign follows from the 
definition of B 1 ' 3 . Since 

n H \ n° H = {t eTi H \ n° H ■. u - tj > -1, v« < j} , 

the part (iii) follows immediately from (ii). □ 

The above proposition provides a decomposition of the boundary into non- 
overlapping boundary elements. To make each boundary clement explicit, we use 
symmetry. Let Q = S4 be the permutation group of four elements. For t G M, H and 
a G Q, the action of a on t is denoted by to~, which means the permutation of the 
elements of t by a. A moment of reflection shows that, for (I, J) C /C, 

B \i\-\J\ _ (J Blja .- {to- : t G Bi tJ , a eg}. 
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Later in the section we will need to consider points on the boundary elements 
that are congruent module H. For (/, J) C JC we further define 

[Bi,j] := {B^j + k: k G Z 4 H } n Ti H = {t + k G H H : t G B/../,k G Z^} . 

Since [-Bj.j] is a subset of £Ih and -Bj,j is a boundary element, we see that [£?/,./] 
consists of exactly those boundary elements that can be obtained from Bi.j by 
congruent modulus H , as confirmed by the following lemma. 

Lemma 3.7. I/t, s G and s = t (mod -ff), t/iert t = s. 

Proof. By Lemma T3.31 if t,s G and t = s (mod H), then s — t G Z|^ and, set 
k := s — t, — f < ki — kj < 1 for alH, j G N4. The last condition means that either 
ki G {0, 1} for all i G N4 or ki G {0, — 1} for all i G N4. The homogenous condition 
ki + k 2 + fc 3 + k 4 — Q then shows that fci = k% = &3 = fei = or s = t. □ 

As an example, we have 
(3 g) B {i},{2,3} = {(t,t - 1, t - 1, 2 - 3t) : i < t < f } , 

fl {1 ,2},{3} = {(l-t,l-t,-t,3t-2): I < t < |} , 
and from the explicit description of i?{i} l {2,3} we deduce 
[5 { i } , {2 , 3} ] = £{i},{2, 3} U (B { i }){ 2,3} + (-1, 1,0,0)) U (B {1} , {2>3} + (-1,0, 1,0)) 

(3.10) = -B{1},{2,3} U -B{2},{1,3} U -B{3},{1,2}- 

Others can be deduced similarly. The last equation indicates that [-B/,j] is a union 
of Bj/ t ji, which we make precise below. 

Let o~ij denote the element in Q that interchanges i and j; then t<7y = t — (t,- — 
tj)eij. For a nonempty set / C N4, define Qi := {a^ : j, j e /}, where we take 
o~ij — Uji and take o~jj as the identity element. It is easy to verify that Qi forms a 
subgroup of Q = S4 of order |/|. 

Lemma 3.8. Let (I, J) G /C. Then 

(3-11) [Bi,j}= (J B ItJ a. 

Proof. For any i, j G /UJ, the definition of Bj j shows that to^ — t = (tj — t^e^ G 
for all t G Bj j. It then follows from toy G that jCTy C [B/j]. 
Consequently, Uo-egjuj ^ t- 8 /^]- 

On the other hand, for any s G [-Br,j] there exists t G -Bj, j such that s — t G 1* A H . 
It follows from Proposition [3761 and Lemma l3~7l that s G Bi l! j 1 for a pair (Ji, Ji) G 
/C. By the definitions of B/j and Sljy , there exist t, s € K such that 

-|<*-l = *j<*j <**=*< f, ie/, j eJ, I^IUJ, 

- I < 5 ! - 1 = Sj < si < Si = s < f , £ e Ji, j G Ji, I & h U Ji. 

Since s — t G Zjj-, the above inequalities imply that Sj — U G {—1, 0, 1} for all i G N4. 
We claim that t = s. Assume otherwise, say s > t. For i G I\, Sj — ti > s — t > so 
that Si~ti = 1; while for i $ Ii, s^— i, > s— 1— t > —1 so that Si—ti G {1, 0}. It then 
follows that X^eN 4 ( Si — - > 0, which poses a contradiction to the homogeneity 
of s — t. Hence we must have s — t. With s = t, it is then easy to see that 
Si — s — Si — ti G {0, —1} for i G /, Sj — s + 1 = Sj — tj G {0, 1} for j G J, 
and s — 1 < si = ti < s for I £ I U J. This shows that I LI J = Ji U Ji and 
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E iG zuJ*i = E tG / lU J 1 s *- Meanwhile, we note that £ iezuJ -ti = *(I J I + l J D ~ l J l 
and £ ie7lUJl s « = s (l J il + I J il) ~ I J il- It then follows that |J| = |Ji| and |7| = |ii|. 
Consequently, s = ta for a a £ and [.Bj.j] C Uo-eSi 7 Bi.jcr. This completes 

the proof of the lemma. □ 

Since A?'- 7 can be obtained from ICq 1 from the action of Q, it follows that 

(3.12) B^ = |J [S/,j]= (J [fl], 0<i,j<i+j<4. 

We also note that [Bjj] n [-Bj^jJ = if (I, J) 7^ (Ii, Ji) for (J, J) £ /C and 
(Ji , Ji ) G /Co , which shows that (|3.12p is a non-overlapping partition. 
If |7| + I J| = 3 or Bj j is an edge, then we have 

B 1 ' 2 = [#{!}, {2,3}] U [B {1] . {2A} ] U [5 { 1}, {3 ,4 } ] U [B { 2},{3,4}], 
B 2 ' 1 = [-B{i )2 },{3}] U [-B{1, 2 },{4}] U [-B{1,3},{4}] U [ B {2,3},{4,}], 

where, recall that Bm/2,3} and ^{1,2}, {3} are given in (|3.9p , 

-B{1},{2,4} = - B {1},{2,3}°'34, ^{1,2},{4} = ^{l^MS} ^! 

(3.14) #{1},{3,4} = -B{1},{2,3}°24, S{1,3},{4} = ■E{1,2} 1 {3}< 7 '230'34) 

^{2}, {3,4} = ^{l,2},{3} cr 120'24, #{2,3},{4} = i?{l,2} : {3} cr 130'34 ■ 

If |/| + |J| = 4, then 

£ 1,3 = [{(|) 3> |> -|)}] ' ^ 2 ' 2 = 2' ~5' ~|})] 

^ a,1 = -?)>]■ 



( 3 - 13 ) 



(3.15) 



3.3. Dodecahedral Fourier partial sum. In order to apply the general result 
on discrete Fourier analysis in the previous section to fee lattice, we choose A = A 
and B = nA with n being a positive integer. Then the matrix 

(2n n n 
n 2n n 
n n 2n 1 

has integer entries. Note that N is now a symmetric matrix so that An = Ajvtr, 
and it is easy to see that A N = An- Recall the definition of H in (|3 . 5|) . Using again 
j = AH(A tr A)~ 1 k £ Zjj, it is easy to see that k £ An becomes j £ H„, where 

i„:= {kei:^e n H } = {k £ H : -An < h - kj < An, 1 < i < j < 4} . 



The finite dimensional space H.n of exponentials in Theorem 12.51 becomes 

Tin := span {0 k : k £ H n } with dimW„ = det(iV) = An 3 . 

Note that the points in H ra are not symmetric under Q, since points on half of 
the boundary are not included. For reasons of symmetry, we further define 

H* := {k £ H : £ U H } = {k £ H : -An < h - k 3 < An, 1 < i < j < 4} . 

For the Fourier expansion (|3.7|) of an 77-periodic function, we define its dodeca- 
hedral partial sum as 

(3.16) S n f(t) := ]T (/, k )0 k (t) = \f f(B)D* (t - g)<fe, 
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(-2n,-2n,2n,2n) 



Figure 3.4. i; = {tel: -An <t t - tj < An}. 



where is the Dirichlet kernel for the dodecahedral partial sum 
(3.17) D*(t):= ]T ^ kt - 



keH* 



Our immediate goal is to find a compact formula for the Dirichlet kernel D% . 
We start with an observation that the index set H* can be partitioned into four 
congruent parts, each within a parallelepiped, as shown in Figures 3.5-3.8. 




(-3n,n,n,n) 



(2n,-2n,-2n,2n) 



(-2n,-2w,2».2n) 




(-n,-„,3»,-») 



Figure 3.5. 



Figure 3.6. 



ir(2) 



(2n, 2n, -In, -2n) 



{2n,2n,-2n,-2n) 




(0,0,0,0) 



(—n, -n, — n,3n) 



Figure 3.7. h! 3) 



(St., -»,-»,-») 




(-2n, 2.1,2?, . -2„) 



Figure 3.8. h! 4) 
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Lemma 3.9. Define M ( n k) := {j G H : < j t - j k < An, I G N 4 } for k G N 4 and 
i„':={kel:i, = fej, Vi,j £ J; and < fc, - k 3 < 4n, Vj £ J, Vi G N 4 \ J} 
/orlcJC N 4 . T/ien 

H;= (J H« and H„ 7 = f| E#). 

Proof. For k G H*, let fcj = minjfci, k%, fc 3 , fc 4 }. Then k G Hl£*\ which implies that 
H; C U 4£N4 Hl j) . Since C H* for each i G N 4 , it follows that H* = [J» eN4 H^. 

If k G n Hl j) , then < h - kj < 0; that is, h = kj. It follows that if 
k G f) jeJ Mti\ then k, = kj, G J, which implies ftjzj^ C H^. Since 
H;( C by definition, we conclude that = f] je J H^. □ 

Theorem 3.10. For n>0, 

4 . 

(3.18) £^(t) = e n+1 (t)-e n (t), wftere 6 n (t) := TT E^IILi, 

*- A Sin7TIi 

Proof. Using the inclusion-exclusion relation of subsets, we have 

D^(t)= E (-i) |J|+1 E e¥k,t - 

0C./CN 4 keH^ 

Fix j G J, using the fact that tj = — J2i^j *» i we have 

^ e* k -t = ^ e ¥ E 1€N4U (fci-fcj)*i = ^ Yl e^( k '- k i )tl . 

keHJ keH^ kGH^(eN 4 \J 

By the definition of and the fact that k G H implies fc, = fcj (mod 4), we obtain 
^ e f' = JJ ^ e *(fci-fcj)*i 

keH^ ZeN 4 \J 0<fc 1 -fc :j <4n 

A:; (mod 4) 

= n e e2 " fc,t; n ^«(*o- 

/eN 4 \J0<fe,<n ;gn 4 \j 

Consequently, we obtain 

(3.19) D%(t)= E (- 1 )' J|+1 II K n (t l )=Y[K n (t j )-l[(K n (tj)-l), 

0CJCN 4 l£N 4 \J i6N 4 jeN 4 

where the second equality is easily verified upon expanding the right hand side 
explicitly Thus, we conclude that 

4 e 2wi(n+l)tj _ J 4 e 27ri (n+l)tj _ e 27ritj 



^)=n -n 



^ -,7rz (n-\-l)tj n — 7zi (n-\-l)t-j 4 „7i"i ni, „ — 7rint-j 

n e e _ e *<n«i _ TT e e J c ^(»+i)t, 

giritj g—Tritj XX gKitj g— 7riij 

4 4 

-p-r sin7r(n + l)tj -pi- smimtj 
sin7rf,- sin7rt,- ' 
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where in the last step we have used the fact that Yij=i e latj = 1, which follows 
from the fact that t\ + t 2 + t 3 + t 4 = 0. This completes the proof. □ 

As an immediate consequence of Theorem 13.101 we conclude that 

|H;j=A?(0) = (n + l) 4 -n 4 . 

The explicit formula of the Dirichlet kernel also allows us to derive an estimate 
for the norm of the partial sum S n f in (|3.16p . Let ||/||oo denote the uniform of 
/ S C(Ojj) and let H^Hoo denote the operator norm of S n : C(£l#) C(Sl ff ). 

Theorem 3.11. There is a constant c independent of f and n such that 

\\S n \\oo <c(logn) 3 . 

Proof. From (|3.16[) . a standard argument shows that the norm is given by 

||Sn||oo = ^ max / \D*(t-s)\ds. 
^ ten H Jn H 

To estimate the integral, we use the first equation of (|3 . 1 9|) and the fact that 



(3.20) 

which leads to 

WSJ 



K n (t)=J2e 

j=0 



2-nijt 



„™( sin7r (" + 1 ) t 



sin irt 



< i max X] / II \ K n(U - si)\ ds 



^2^ x ^2 i n 

EUh 0C,7CN 1 ,/Off ZeN 4 \J 



sin7r(n + l)(ti - si) 



sin7r(t; - s t ) 



ds. 



Since J ^ 0, the above product contains at most three terms, and those that 
contain product of three terms dominate other integrals. Consequently, enlarging 
the domains of the integration and then using the periodicity of the trigonometric 
function, we conclude that 



\Sn\ 



< C 



n 



■i,i]» J= i 



sin7r(n + 



sin mij 



du 



sin7r(n + l)Uj 



sin 7ruo 



duj < c(logn) 3 



where the last step follows from the usual estimate of the integral involved. 



□ 



We expect that the estimate is sharp, that is, ||SVi||oo > c(logn) 3 . To prove such 
a result would require a lower bound estimate of the integral of \D^(t — s)| at one 
point in S7, likely at s = 0. However, this does not look to be an easy task as there 
is a sum of four terms of the same type. 
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3.4. Discrete Fourier analysis on the rhombic dodecahedron. Using the 
set-up in the previous subsection, Theorem 12.51 in the homogeneous coordinates 
becomes the following proposition 



Proposition 3.12. For n > 0, define 

1 

in 3 



J 61 



Then 



</, 9) = </, 9)r 



f, 9 G Ur. 



The point set H n , hence the inner product (•,•)«, is not symmetric on in the 
sense that it contains only part of the points on the boundary. Using the periodicity, 
however, we can show that the inner product (•, •)„ is equivalent to a symmetric 
discrete inner product based on Ml. To proceed, define 



and, recall 
(3.21) 



ii'.j 



K-= {jGH:iefi^} 
for < i, j < i + j < 4 define 
:={kGH: T77- />""}• Kf,0 ■= {k G 



4n fc -°0 J 



Recall that .B^ is a boundary element of £Ih, so that describes those points 
j in H„ such that £ are in of <9ft ff . Furthermore, = {B™ a : a G Q}. 
Using Proposition GOB it is easy to see that Wtf fllj- 1 = if i ^ k, j ^ Z, 

(J H^'=H;\H° and |J ] 

0<ij'<i+j<4 0<i,j<i+j<4 

Lemma 3.13. For n > 1, < i,j < i + j < 4, 

4! 



f,3 



H„\H°. 



(n-1) 4 



;i ;iu - ; - j)\ 



-An-l) 4 -^. 



Proof. The first equation follows from |H°| = |H* 
tion of B 1 ^ in Subsection 3.2 shows that B 1 ^ has 



4! 



- (n — 1) . The descrip- 
segments and each has 
□ 



(n — l) 4 1 J points, which proves the second equation. 
Definition 3.14. For n > define the symmetric discrete inner product 



j6H* 



ere c ; = 



1 i/j G H°, and Cj = ^ttj */j G M^f ; more explicitly 



{n) 



1 

2 ' 
_L 

3 ' 
_L 

4 ' 
1 



(ri 4 — [n — l) 4 points in the interior), 

(12(n — l) 2 points on the faces), 

j E JHL£* U H 2 ' 1 , (2 x 12(n - 1) points on the edges), 

j G H*' 3 U H^ 1 , (2x4 points on the vertices), 
rr2,2 



G 1H1„, 

j g ebi, 

[Tl,2 



^ ... j G (6 points on the vertices). 

It is easy to verify that J^jeH* c j™^ = ^" 3 ' so ^' ^» = ■"" ^ e P rove tne 
following result. 
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Theorem 3.15. For n > 0, 

(/, g) = (/, g)n = (f, g)* n , f,ge n n . 

Proof. For each j G H^q , we define 

(3.22) 5j={kei;:i=l (modi/)}. 

It follows immediately from (|3~12|) that Wtf = UkeH^ % HSHJ, [-B 7 " 7 ] is the 
union of (' 7 '|j|' J ') components oiBj jcr. Consequently, it follows that |6j| = ( fe j~ ( ) 
for j G H^q. Let / be an ii-periodic function. Then 

jGH*\H° 0<i,fc<i+fc<4 V i / jgg*.* 



Since |iSj| = ( ! ^ fe ) for j G M^ k , using the invariance of /, we then conclude that 



c 



jeH*\H° 0<i,fc<i+fc<4 V i ) jgH i, '" a jeH„\H° 

Since cj = 1 if j G H° , the proof is completed. □ 

The discrete inner product is closely related to cubature formula, since Theorem 
13.151 shows that the integral of / G H. n agrees with the discrete sum over H* . 
In fact, more is true. Let us define by T n the space of generalized trigonometric 
polynomials, 

T„:=span{0 k :keHj;}. 
Theorem 3.16. For n > 0, the cubature formula 

1 f ,. 1 \- („). 



/( t )dt = > c\ n >f(±) 



is exact for all f G 72 n _ 



Proof. If k,j G Bin, then the definition of implies immediately that k - j e 
H| n _ 1 . Suppose now j G and we may assume that ji > j% > jz > ji- 

There exists k G H* such that k j G H* . Indeed, if j\ - h < 4n - 4, we can 
take k = (n, n, n, — 3n). If j2 — J4 < 4n — 4, we can take A; = (3n, — n, — n, — n). 
Finally, if both ji — j'3 > 4n and j'2 — J4 > 4n then it follows from the definition 
of H^.-l that ji — j 2 < 4n — 4 and j'3 — j'4 < 4n — 4. In this case we can take 
k = (2n, 2n, — 2n, —2n). Consequently, this shows that 

HS B _ 1 = {l:l = k-j,k ) j6l B }. 

Thus if 4>i £ ^2n-i, then j G H2„_ 1 and there exist k, 1 G H n such that (j>$ — <^k</>i- 
Consequently, the stated result follows from Theorem 13. 151 □ 
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3.5. Interpolation on the rhombic dodecahedron. For the rhombic dodeca- 
hedron, Theorem 12.61 on interpolation becomes the following: 



Proposition 3.17. For n > 0, define 

I n f(t) := f(±)*n(t - jL), where $„(t) = JL £ 0fc ( t ), 



for f £ C(CIh )■ TTien !„/ G 7i„ and 

= Vj G H„. 

Again there is a lack of symmetry in the sense that X n uses only points in H n , 
which contains only about half of the boundary points. We are more interested in 
another interpolation operator given below, defined over all points in H* . Although 
it does not interpolate at all points in H* , its symmetric form can be used to derive 
results on the tetrahedron in the next section. Recall 5k defined in (|3.22[) . 

Theorem 3.18. For n>0 and f G C(TLh), define 

Km ■■= E /(£)<j,»(t), 



wht 



*j, n (t) = *;(t - and S* (t) = JL J2 4"Vk(t). 



keH* 



Then T*f £ T n and it satisfies 
(3.23) T*f(±) = 



f ( in J ' 



j e 



Furthermore, $^(t) is a real function and it satisfies 



An 3 



(3.24) 



Proof. By definition, 



(£.^(t)+^_ 1 (t))--E 



1 sin(n — l)7ri„ 



j/=i 



sin 7rf „ 



cosri7r(2tj + i,/) 



j'=i 



1 4 1 

- >J cos 2irntj — - cos 2'im(t fjt + £„) 



3=1 



l<fi<v<A 



It] 



Since tiles Rjj-, there exist m, 1 G 1 A H such that G and k — j = m + 4nl. 
Thus, by Theorem 13. 151 

=i E = ^E c ^(ir) 

iGH* ieH* 
= (0m) 4>o)*n = (<f>m,<M = <$m,0- 
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Equivalently we can write the above equation as 



(3.25) *J,»(£) = fak>4l>* 



1, k = j+4nl,leZk 
0, otherwise, 



the values of that 



from which (|3.23[) follows. 

To derive the compact formula for iy n we start with the obvious fact that <9H* = 
H* \H*_i, so that 

E <4 B) *k(t) = 5 E Mt)- E (l-4 n) )^(t) 

kedH* keHJ\B£_i kGOH* 

= i(Df(t)- J D« 1 (t))- e (i-4 B) )^(t). 

kedH* 

Since H* = H° Ucffl* , we then derive from the decomposition of cffl* into and 

l -( D X{t) + D*_ x {t))- l ~ E Mt) 

fceH^UH 2,1 

-\ E Mt)-| E m*) 

fceH^uH 3,1 fceH 2,2 

Let us define H^ J := {k e H : i e for /,JcN 4 and also define [H£ J ] := 

{k e I : i 6 [-B/,j]}. It follows from ([3~T2]) . (|3~2Tj) and Lemma iU that 

h# = U [<' J ] and t H « J ] = U H « v 

In particular, by (|3 . 13[) . we have 

H,V 2 = [MW'W] U [BUi 1 ^ 2 ' 4 }] U [H^ 1 ^'^ 3 ' 4 ^] U p4 2 M 3 .*>], 
H 2 ' 1 = ppM 3 )] U {MV&'W} U [Hp^ 4 *] U [H{2' 3 >-{ 4 >]. 

By (BD , [ui 1}A2 ' 3} ] = ui 1}A2 > 3} U H£ 2} ' {1 ' 3} UH< 3Ml ' 2} . Furthermore, it fol- 
lows from j3H that [H< 1} ' {2 ' 4} ] = [H< 1} ' {2 ' 3} ]a34, [H< 1} ' {3 ' 4} ] = [H UM2 ' 3} ]a 24 , 
[Hl4 2 ^'^ 3 ' 4 ^] = [IHlii 1 ^'^ 2:3 ^]o'i20'24- Using the explicit formulas in (|3.9[) . it is easy to 
see that 

H {i},{2,3} = {(j + 2n, j - 2n, j - 2n, 2n - 3j) : 1 < j < n - 1} . 
Consequently, using t\ + 1% + + ti = 0, it follows readily that 

E &oo= E ^(t)+ E ^(t)+ e m*) 

ke[ni 1} - {2 ' 3} ] kGH* 11 ' 12 ' 31 keH< 2 >< fl ' 3 > kGHi 3 >' {1 ' 2 > 

n-1 

_ y ^ e ~2Trijt 4 ^ e 2n7ri(ti+t 4 ) ^ ^nmfo+ti) _|_ e 2ri7ri(t3+t 4 ) N 



sin(n — l)7rt4 



^ e n7ri(2ti+t 4 ) _|_ e ri7ri(2t 2 +i 4 ) ^ e ri7ri(2i 3 +t 4 )^ 



sin 7rf4 

where in the last step the sum is evaluated using (|3.20p . The explicit formulas for 
other components of X^keH 1 - 2 ^ k follow from the above expression by permuting 
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the variables. In a similar manner, we have [H^ 1 ' 2 ^ 3 ^] = H^ 1,2 ^ 3 ^ U H^ 1 ' 3 -^ 2 ^ U 
HF } ' {1} , and, using &® again, [Hp } ' {4} ] = [Hp } ' {3} ]a 3 4, [Hp } ' {4} ] = 
[H^' 2},{3} ]o-230-34, [Il£ 2 ' 3} ' {4} ] = [H^'^'^Vis^- Moreover, we also have 
H {i,2},{ 3 } = {(2n _ i; 2n _ i; _ 2n _ jf _ 2n + 3j) :l<j< n -l}. 

Thus, using t\ + ti + ts + <4 = 0, we can deduce as before that 

J2 E E e k (t) 

kepi 1 ' 21 -* 31 ] keHi 1 - 21 ' 131 keiii, 1 - 3 ^ 21 keHi 2 ' 3 >- {1 > 



sin(n — 1)77^4 
sin Ttt^ 



^ e -n7ri(2ti+t 4 ) _|_ g — n7ri(2t 2 +i4) _|_ g-n-n-i(2t3+t 4 )^ 



from which the explicit formulas of other components of X)keH 2,1 ^ k f°U ow from 
permuting the variables. 

Putting the sums over H*' 2 and H 2 ' 1 together, we obtain 



E, , . „ v-^ sinfn — l)7ri„ v-^ ,„ x 

<f> k (t =2> — i '— -> cosmr (2t 1 + t v ). 
^— ' sin7ri,. 



3=1 

Using (|3.15[) . it is easy to see that H^ 3 = {(n, n, n, — 3n)<7 : cr e 5}, HI 2 ' 2 = 
{(2n, 2n, — 2n, — 2n)cr : <t e 5} and H 3 ' 1 = {(3n, — n, — n, — n)a : a £ G}, from 
which it follows that 

4 4 
^ fe (t) = ( e2mt3 + e~ 2mt i ) = 2 ^ cos 2^ . 

Furthermore, using t\ + <2 + £3 + ^4 = 0, it is easy to see that 

Mt)= e 2m[tM = J2 coB27rn(t M +t w ). 

kGH 2,2 l</i<i/<4 1</j<i/<4 

Putting these terms into (|3.26[) completes the proof. □ 

The compact formula of the interpolation function allows us to estimate the 
operator norm of I*, which is usually referred to as the Lebesgue constant. 

Theorem 3.19. Let ||2£||oo denote the operator norm of I* : C(Clu) 1— » C(Qh). 
Then there is a constant c, independent of n, such that 

II^JU <c(logn) 3 . 

Proof. A standard procedure shows that 

Halloo = max ]T |$;(t-^)|. 
ten " ken* 

Using the compact formula of $* in Theorem 13.181 it is easy to see that it suffices 
to prove that 

-Lmax E \ D n(t-M <c(logn) 3 , n>0. 

m te ° H k6Hl 
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Furthermore, as in the proof of Theorem EUil the formula of D% in ([3l9]) shows 
that our main task is to establish the estimate 



/{i,2,3} := T~~3 max £ \K n {t x - ^)K n (t 2 - &)K n (t. 

4 " te °"kGH* 



., - ha.) 

J An' 



< c(logn)" 



and three other similar estimates 1/1,2,4}) ^{1.3,4} an d ^{2.3,4} , respectively. Enlarg- 
ing the domain H* to {k S 1 A H : — 4n < fcj < 4n, fe, = (mod 4), 1 < i < 3}, we 
see that 



j n n <l 

/ { i,2,3}<^ max V E E l^fe-^-^-l) 



te[-i,i] 



fci— — n ko — — n. Ai3 — — n 



1 / 1 

< — max — > 

4 teT— 1,11 \ n 



fc=-r 



sin(n+l)7r(t-f) 



sin7r(i - f) 



< c(logn) 3 



where the last step follows from the standard estimate of one variable (cf. [131 Vol. 
II, p. 19]). □ 

Again we expect that the estimate is sharp, that is, ||/*|| > c(logn) 3 ; and the 
problem is again that there is a sum of four terms of the same type. 



4. Discrete Fourier analysis on the Tetrahedron 

Considering functions invariant under the isometrics of the fee lattice, the dis- 
crete Fourier analysis on the dodecahedron in the previous section can be carried 
over to the analysis on the tetrahedron. 

4.1. Generalized sine and cosine functions. The fee lattice is the root lattice 
of the reflection group A3 Chapt. 4]. Under the homogeneous coordinates, the 
group A3 is generated by the reflections joy : 1 < i < j < 4}, where cry is the 
reflection defined by toy = t — 2 ^j t,e ' e J ^ with e^- = — ej as before. Thus 
.A3 is the permutation group in the previous section. Denote the identity element 
in A3 by 1. It is easy to see that we have 

= l, OijOjkVij = &ik, i,j,k e N4. 

For o G Q = A3, let \a\ denote the number of inversions in a. The group Q 
is naturally divided into two parts, Q + := {o G Q : |o| = (mod 2)} of elements 
with even inversions, and Q~ := {o G Q : |c| = 1 (mod 2)} of elements with odd 
inversions. Writing it out explicitly, we have 

G + = {l, cri2(Ti3, 0-13CT12, cr 12 CTi4, (Ti4(Ti2, fTi3(Ti4, 

C14CT13, O23CT24, 024 23, O12O34, O13O24, 0-14(723} 
Q~ = {012, O13, O14, O23, 2 4, O34, O12O13O14, O12O14O13, 

013012014, O13CT14O12, 0-14O12CT13, O14O13O12}. 

The action of o G Q on the function / : IR^ i — > 1R is defined by o/(t) := /(to). A 
function / in homogeneous coordinates is called invariant under Q if of = / for 
all o G Q, and it is called anti- invariant under Q if af — p{o)f with pier) = 1 if 
o G Q + and p(o) = -1 if o G Q~ . 

The following proposition follows immediately from the definition. 



22 



HUIYUAN LI AND YUAN XU 



Proposition 4.1. Define two operator V + and V acting on /(t) by 



(4.1) 



O ± 



<?eg- 



/(tcr) 



T/ien i/ie operators V + and V are projections from the class of H -periodic func- 
tions onto the class of invariant, and respectively anti- invariant functions. 



Applying the operators V ± to 0k (t) 



tk-t 



gives basic invariant and anti- 



invariant functions, which we denote by TCk and TSk, respectively, as they are 
analogues of cosine and sine functions. We formerly define them as follows. 



Definition 4.2. For k € H define 

TC k (t) := 7>+0 k (t) = ^ 



<^k(tcr) + ^2 M t<7 ^ 

<ree+ aeg- 



TS k (t) := -7>"0 k (t) = ~ 



creG + 



0k (to-) 



aeg- 



0k (tcr) 



and call them generalized cosine and generalized sine, respectively. 

Evidently TCk is invariant and TSk is anti- invariant. The rhombic dodecahedron 
is invariant under our group Q of order 24; its fundamental domain is a tetrahedron. 
For invariant functions, we can make use of symmetry to translate results on the 
rhombic dodecahedron to one of its 24 tetrahedrons. We shall choose our reference 
tetrahedron as 

(4.2) A:= {ieK 3 :0<i 3 ±i 2 ,i2±J;i<l}. 

In the homogeneous coordinates, by (|3.2[) , this tetrahedron becomes 

A H := {t E R% : < t x - t 2 , h - t 3 , t 3 - U, t x - U < 1} . 

See Figure 4.1 below, in which coordinates of the corners are given in both ffi 3 
coordinates and homogeneous coordinates in R 4 . 




(0,0,0.0) 
=(0, 0,0) 



Figure 4.1. Reference tetrahedron. 



When TCk are restricted to the tetrahedron Ah, we only need to consider a 
subset of k £ H. In fact, it is easy to see that 

(4.3) TCkcr (t) = TCk(to) = TCk(t) 
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for t G Ah and a G Q. Thus, we can restrict k to the index set 

A := {k e H : ki > k 2 > k 3 > fc 4 } 

when we consider TCk- As for TSk, it is easy to see that TSk CT (t) = TSk(ter) = 
TSk(t) for a G Q + and TSko-(t) = TSk(tcr) = — TSk(t) for a G G~ ■ In particular, 
TSk(t) = whenever two or more components of k are equal. Thus TSk are defined 
only for k G A , where 

A° := {k G H : k x > k 2 > k 3 > k 4 } , 

which is the set of the interior points of A. Since k G H implies that ki G Z and 
ki + k 2 + k 3 + ki = 0, the points in A lies in a three dimensional wedge. To describe 
the points on the boundary of A, we further define 

A-^ := {(k G H : k\ = k 2 > k 3 > k^orki > k 2 — k 3 > k^orki > k 2 > k 3 = k^} , 

A e ' 2 := {(k, k, k, -3k), (3k, -k, -k, -k) : k > 0} , 

A"- 1 := {(2k, 2k, -2k, -2k) : k > 0} , A" := {(0, 0, 0, 0)}. 

Then evidently A \ A° = A' U A 6 ' 1 U A e < 2 U A' u . 

Let kQ denote the orbit of k under Q, that is, kQ :— {ka : a G G}- Then, for 
k, j G A, kQ fl jQ = whenever k ^ j. Furthermore, it is easy to see that 

'24, k G A°, 
12, keA-f, 



(4.4) 



TCk(t) 



1 



|k<?| = < 



G. 
4, 
1, 



k G A 6,1 , 
keA e > 2 , 
k = G A v . 



We define an inner product on Ah by 
1 



(/,<?> 



IA 



H 



A, 



/(t) fl (t)dt = 12 / f(t)g(t)dtidt 2 dt 3 . 



A/ 



If fg is invariant under Q, then it follows immediately that (f,g) = (f,g)A H - 
Furthermore, the generalized cosine and sine functions are orthogonal with respect 
to this inner product. 

Proposition 4.3. For k, j G A, 



(TC k ,TCj)A H = k '' 



= Sui < 



(4-5) m 
for k,jGA°, 

(4.6) (TS k ,' 

Proof. Both of these relations follow from the identity (f,g) = (f,g)A H f° r invari- 
ant functions. For (|4.5[) . the invariance is evident and we only have to use the 
orthogonality of 0k in Proposition 13.11 and (|4.4[) . For (|4.6[) , we use the fact that 
TSk(t)TSj(t) is invariant under Q and the orthogonality of </>k on f2#. □ 



1, 


k = 0, 


1 

4' 


k G A e > 2 


1 

6' 


k G A 6 ' 1 


1 

12 ' 


k G A f 


1 

> 24 ' 


k G A°; 


1 k ' J 
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We can derive more explicit formulas for the generalized cosine and sine functions 
by making use of the homogeneous coordinates. For example, we have 

4>k(tl, t 2 , t 3 , ti) — (f>k{t2,tl,ts, t 4 ) — <pk(ti,t2, t 4 , t 3 ) + 0k(i2, ti, t 4 , t 3 ) 
_ ^ e ^(fciti+fc 2 t 2 ) _ e ^(fcit2+*2ti)j ^ e ^-(k 3 t 3 +k 4 t 4 ) _ e ^-{k 3 t 4 +k 4 t 3 )j 

= 2ie^ (fcl+fc2)( * 1+ * 2) sin ^(^-^(ti-^) . 2ie^ (k3+ki)(t3+t ^ sin ^-Mfe-U) 

= -4 e T[C £ l+ fc 2)(tl+t2) + (fe 3 +fc4)(t3+*4)] sm ^{kl-k 2 )(tl-t 2 ) ^ ir(k 3 -k 4 )(t 3 -ti) 

= 0k(^3, *4, tl, £2) — 0k(t4,*3)*l)*2) - 0k (*3, *4, *2, *l) + 0k(*4, *3> *2, 

Similarly or by permuting the variables, we also have, 

0k(*l,*3,*4,*2) - 0k(*3,*l)*4,*2) - 0k(*l,*3,*2,*4) + 0k(*3>*l>*2,*4) 
= <f>k(t 4 ,t 2 ,h,t 3 ) - 0k(*4,*2,*3,*l) - 0k(*2,*4,*l,*3) + 0k(*2, *4, *3, *l) 

= -4 e T i [( fe l+ fc 2)(tl+t3) + (*:3+fc4)(i2+i4)] g m ^(fcl-fc2)(tl-«3) gm (fcg -fc 4 ) (t 4 - 1 2 ) 

and 

0k(*l)*4) *2j ^3) - 0k (*4, il, *2, £3) ~ 0k(*l>*4j *3) £2) + 0k(*4> *1> *3> *2) 
= 0k(*2)*3) *1> £4) - 0k(*2, ^3, £4, £l) - 0k(^3, *2j tl> £4) + 0k(*3, h, t 4 , h) 

= -4 e ¥[( fe l+ fe 2)( t l+ t 4) + (' £ 3 + fe 4)(*2+*3)] s i n "-(fcl-fc2)(tl-t4) gm (k 3 -kj) (t 2 - 1 3 ) ^ 

Consequently, using the homogeneous relations t\ + t 2 + t 3 + t 4 = and k 4 + k 2 + 

1 2 ) sinf (k 3 - k 4 )(t 3 - t 4 ) 

1 3 ) sinf (fc 3 - fc 4 )(*4 - h) 
t 4 )sin-(fc 3 - k 4 )(t 2 - t 3 ). 

1 2 ) cosf(fc 3 -k 4 ){t 3 -U) 

1 3 ) cos|(fc 3 -k 4 )(t 4 -t 2 ) 

1 4 ) cosf(fc 3 -k 4 )(t 2 -t 3 ). 

Permuting variables ti,t2,t 3 , t 4 lead to other representations of TSk and TCk- 

4.2. Discrete inner product on the tetrahedron. Using the fact that TCfc and 
TSfe are invariant and anti-invariant under Q and the orthogonality of fa with re- 
spect to the symmetric inner product (•, ■)* , we can deduce a discrete orthogonality 
for the generalized cosine and sine functions. For this purpose, we define 

(4.9) A„ := H; n A = {k G H : k 4 < k 3 < k 2 < fci < k 4 + An} . 

The point set with n = 4 and the region is depicted in Figure 4.1. 



k 3 + k 4 = 0, we can then deduce that 

(4.7) TSk(t) =i e *(fci+fc2)(ti+t2) sin f (fci - k 2 ){h - 

+ i e ^( fel+fe2 )( tl+ t 3 ) sin _ k2 ^ h _ 

+ i e¥(fcl+fe2 ) (tl+u) gin |( fci _ k2 ^ h _ 
In a similar way, we obtain 

(4.8) TCk(t) =l e ¥(ki+k 2 )(t 1+ t 2 ) cos - ( fcl _ Wl _ 

+ i e ¥(fei+fc2)(ti+t 3 ) cos |( fcl _fc 2 )( tl - 
+ i e¥(kl+k2)(tl+ti) cos | (fci _ k2 ^ h _ 
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(2n,2ra,-2ra,-2rc) 
= (0.0,4*.) 




(O.O.I). 0) 
= (0.0,0) 



Figure 4.2. A„ with n = A 



By the definition of H* , the set : k e A„} contains points inside Ah- Wc 
will also need notation for points on the boundary of Ah, which are defined as 

A° := H° n A°, Ai := (A' n H°) U (A° n H^ 1 ) , 

A^:= (A^niXA^nHy), 

A^ 2 := (A e < 2 n H°) U (A^ n H^ 2 ) U (A* n H 2,1 ) , 

A^ := {0} U (A 6,1 n H 2,2 ) U (A e < 2 n H^ 3 ) U (A e < 2 n e 3 - 1 ) , 

corresponding to the set of the interior points, the set of the points on faces, the set 
of points on two type of edges, and the set of vertices, respectively. More precisely, 
these sets are given explicitly by 

A° = {k e H : fc 4 < k 3 < k 2 < fci < k 4 + An} , 

A{ = {k £ i : fc 4 < k 3 < k 2 < fci = ki + An or fc 4 < fc 3 < k 2 = k\ < fc 4 + An 

or fc 4 < k 3 = k 2 < k\ < ki + An or fc 4 = k 3 < k 2 < k\ < fc 4 + 4n}, 
A^' 1 = {(2k, 2k, -2k, -2k), (2k + n, n - 2k, n - 2k, 2k - 3n) : < k < n} , 
A^ 2 = {(k,k,k, —3k), (3k, —k,—k,—k),(n+k,n + k,n — 3k,k — 3n), 
(3n — k,3k — n, —n — k, —n — k) : < k < n}, 
A£ = {(0,0,0,0), (2n,2n-2n,-2n), (3n, -n, -n, -n), (n, n, n, -3n)} . 

We denote by TC„ and TS n the spaces of the trigonometric polynomials 

rC n := span {TC k : k 6 A„} , TS n := span {TS k : k e A° } , 

respectively. We define a discrete inner product (•, -)A,n by 
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where 



(4.11) 



(n) 



24, jeA°, 

12, jeA/, 

6, j G A* 1 , 

4, j G A-- 2 , 

1, j G A* 



Theorem 4.4. For /g G TC^n-i, 

(4-12) (/,<?>A„ = (/*<?> A,«. 

Moreover, the following cubature formula is exact for all f 6 7"C2 n -i, 



(4.13) 

In particular, 
(4.14) 



(TC k) TC 



•j/A,n - (n) 
Ai, 



, k,j G A r , 



Proof. We shall deduce the result from Theorem l3.16l Let / be a function invariant 
under Q. Recall the coefficients in the symmetric inner product defined in 
Definition 13.141 Taking into consideration of the orbits of the points in various 
regions, we obtain 

E4 n V(i) = E/(ak) = 24 £ f(D + 12 E /(afc) 

jeH° jgh° jeA°nH° jeA^nH" 

+6 E /(i)+ 4 E /(£)+/(<>). 

jeA^nB" jeA». 2 nH° 

and, using the values of Cj , 

E «, (B) /(£) = m E 4 n) /dk) + i2 E c l n) /(i) 

jgh*\h° jeA^nHi- 1 jeA'nHi- 1 

+ 12 E 4 ,l) /(i)+ 6 E c i n) /(i) 

jeA^nBi 2 ,' 2 



je(A/nHj 1 ' 2 )u(A/nH 2 ' 1 ) 

E 

je(A'=- 2 nHi' 3 )u(A=. 2 nHj 1 ) 



12 E /(ir) + 6 E /(ii) + 4 E 



/ ( 4n 



+ E /(£) + E 



j 6 ( A/ nHft )U(A^ nH„' ) 



/ ( 4 n ) • 



je(A=. 2 nH^ 3 )u(A=. 2 nH^ 1 ) 
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Adding these two expressions together and use (I4.10j) . we conclude that 

(4-i5) E ^nl) = E f(D + E c i n) /(i)- 24 E/(i) 
+ 12 E /(£) + 6 E + 4 E /(£) + E /(£) 

jeAl jeA^' 1 jeA^' 2 JeA- 

= E A j n) /(i). 

jeA„ 

Replacing / by fg, we have proved that (/, g)* L — (/, g) A,n whenever /<? is invari- 
ant. Hence, (|4. follows from Theorem l3.l51 Furthermore, since Jq h f(t)dt = 
i ^ i J Ah f(t)dt for all invariant /, (|4. 13|) follows from Theorem 13. 161 

Furthermore, replacing / by TC k TCj in (|4.15p . we derive by (|4.3| that 

(TC k ,TC j ) A ,„ = * ^ c, (n) TC k (^)TQ(X) 
leH* 

1 1 x Wv^ , , i „ ; i ; 1 1 



4/l :, , 4 E ^ E fc(fc)TCj(£) = ^2lE <T E ^(fe)TQ(fe) 
igh* o-eg leH* o-ee 

IGH* o-GS 

Using (|3.25p and abbreviating k = j mod 4Zj| as k = j, we further deduce that 



1 



24 



TC k , TCj) A) n = 7^7 W G 5 : = k} = -i£ {a e £ : k<r = k} 
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where the last equality follows from a direct counting. □ 

The proof of the above theorem also applies to /, g £ TS*™ since is invariant 
if both / and g are anti- invariant. Moreover, /<? £ Ti n ~\ if /j<? G %i- Notice also 
that TSk(^) = when j G A„ \ A° , we deduce the following result. 

Theorem 4.5. Let the discrete inner product (v)a°,ti ^ e defined by 

g ^ , ; 

(/,3)a°,« = 7^3 E fiidai-k)- 
jeA= 

T/ien 

(f,g)A°,n = </,5 , }a„, /,geT5„. 

4.3. Interpolation on the tetrahedron. We can deduce results on interpolation 
on the tetrahedron by making use of the orthogonality of generalized trigonometric 
functions with respect to the discrete inner product, as shown in our first result 
below. Recall the operator V defined in (|4.ip . 



Theorem 4.6. For n > 0, and f e C(A#) ; define 

C n f(t) := £ /(£)/?, n (t), $ in (t) := ^ E TS k (t)TS k (i). 
jeA° k eA° 

TTien £ n / is the unique function in TS n that satisfies 
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Furthermore, the fundamental interpolation function £? n is real and satisfies 
where means that the operator V~ is acting on the variable t and 0„ is defined 

in (E3HD. 

Proof. By (jH))) . (TSj, TS k ) A °,n = ^^j,k, which shows that ^? n (^) = <5j,k and 
verifies the interpolation condition. It follows from the definition of TS k that 

$.n(t) = ^t-^r E w*)wE- 

keH° 

Furthermore, we can replace the summation over k S H° by the summation over 
kei„ since TS k (^.) = TSj(^) = whenever k 6 A„ \ A°, and TS k (t) = 
whenever at least two components of k are equal. Consequently, we conclude that 

where is the Dirichlet kernel for the rhombic dodecahedral Fourier partial sum 
defined in (|3.17|) . Recall that Q is a permutation group and \a\ denote the number 
of inversions in a £ Q. Let / be an invariant function under Q. Then 

v t -v-f(t - s) = e E(- 1 ) |,r|+|r| /( tff - sr ) 
= 1 4EE(- 1 ) H+|r, /( toT - 1 -) 



^(-l)l CT l/(ta-s) = 7' t -/(t-s) 



where in the third equal sign we have used the fact that err + |r| = \o~\, which can 
be easily verified. Setting / = completes the proof. □ 

The function C n f interpolates at the interior points of A„. We can also derive 
an analog result for interpolation on A„ by using the same approach. However, it is 
more illustrating to derive it from the interpolation on the rhombic dodecahedron, 
which we carry out below. 

Theorem 4.7. For n > and f e C(A H ) define 

_ A (n) 

C*J(t) := E ■= inT E AL' l) TC k (t)TC k (i). 

j£A„ keA„ 

Then L* n f is the unique function in TC n that satisfies 

= /(£), J G A„. 
Furthermore, the fundamental interpolation function £^ n is given by 

*f n (t) - Af ) P+£ j ,„(t). 
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Proof. It follows from (|4.14|) that if n {j^) = #k,j for k,j 6 A„, which verifies the 
interpolation condition. Furthermore, in the proof of Theorem 14. 4[ we established 
that J2jeA n A j #(4^) = Z^eHj ^sK^) for function </ invariant under 0- Apply- 
ing this relation to .g(^) = TCk(t)TCk(^), we obtain 
\(«) 



=4^3 E 4" J TC k (t)TC k (^) 

ken; 



keH* 



4l! - 



Using the fact that Q is a permutation group, it is easy to see that 

v+v+f(t - s) = r+f(t - S ) 

for an invariant function /. Consequently, 

#»(*) = ^KK(* - i) = Aj"¥+£j,„(t). 
The proof is completed. □ 



Recall the explicit formula of given in Theorem 13. 18( £j n enjoys a compact 
formula. 

Let ||Ai|| and ||£^|| denote the operator norms of C n and £*, respectively, both 
as operators from C(Ajj) * C(Ah)- From Theorems 14.61 and 14.71 an immediate 
application of Theorem 13 . 191 yields the following theorem. 

Theorem 4.8. There is a constant c independent of n, such that 

\\£ n \\ < c(logn) 3 and ||£*|| < c(logn) 3 . 

4.4. Interpolation on the regular tetrahedron. The results in the above are 
developed in homogeneous coordinates. Here we indicate how they can be recast 
into the usual coordinates on the regular tetrahedron A* defined by 

A* := {x e R 3 : < x 3 < x 2 < x x < 1} 

as depicted in the Figure 4.3 below. 

The change of variable from t to x £ M 3 is given in (|3.2p . When we transform 
the formulas from the homogeneous coordinates to the regular coordinates, we also 
need to transform the indices from j € Z^ to Z 3 by using 

(4-16) k=-A«U tr j <=► jfc 2 = i(j 2 -j 4 ) 

Under this change variables, it is easy to see that the point set A„ becomes {k € 
Z 3 : < k 3 < k 2 < ki < n}. For example, the cubature formula in (|4. 1 3[1 becomes 
the following: 

Theorem 4.9. For n > 0, the cubature formula 

6 [ J(x 1 ,x 2 ,x 3 )dx 1 dx 2 dx 3 ^^ E ^W^'^t) 

JA * 0<fc 3 <fc 2 <fci<n 

is exact for all f 6 TC 2n -i, where A^." fe2 fc3 = Aj with j given by (|4. 16[) . 
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Figure 4.3. Regular tetrahedron 



We also note that the Dirichlet kernel in (|3. 18|) can be recast into x coordinates 
by (|3.2p straightforwardly, so is the fundamental interpolation function €j in (t) given 
in Theorem 13.181 Consequently, the Lagrange interpolation function in Theorem 
14.71 becomes 

0<fe 3 <fc 2 <fci<ri 

where £ kn satisfies the compact formula 

<» = Al n) P + 4n(t), where £j, n (t) = *;(t - £) 

with $* given in (13"M . A<f } = Aj n) with j as in (|4.16|) . In the above formula we 
apply V + to the compact formula of $* first and then use (|3.2p to change from t 
to x. 

Acknowledgment. The authors thank an anonymous referee for his careful reading 
and invaluable comments. 
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